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It is demonstrated that the fragmentation functions at large momentum fraction play a key role 
in hard hadron production from relativistic proton-proton collisions. We find that this region of 
the fragmentation functions is not strongly constrained by the electron-positron data. This freedom 
can be used (together with the transverse momentum distribution of partons) to reproduce hard 
pion-to-proton ratio data in relativistic proton-proton collisions. 
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Hadron spectra at large transverse momentum (pt), 
where perturbative Quantum Chromodynamics (pQCD) 
has good predictive power, are very important for our 
understanding of the physics at the Relativistic Heavy 
Ion Collider (RHIC) and at the planned Large Hadron 
Collider (LHC). Various experiments at these colliders 
focus on hard (high px) hadron spectra. Jet "tomogra- 
phy" (the study of the strongly-interacting medium via 
energy loss of hard partons) has been proposed to detect 
the quark-gluon plasma (QGP), using high-p^ hadron 
production 



Suppression of total charged hadron production in 
Au + Au collisions relative to a nucleon-nucleon ref- 
erence is reported at RHIC However, for the 
proton-to-pion (p/ir) ratio the PHENIX experiment re- 
ports an anomalous enhancement in Au + Au collisions at 
\/s = 130 GeV ||. An explanation for the p/ir enhance- 
ment was proposed recently, combining pQCD with soft 
physics and jet quenching Q. It should be kept in mind 
in this context that, while pQCD is quite successful for 
total charged hadron (h + + h~) and pion production at 
large px in pp collisions, proton production in pp is not 
well understood using the language of pQCD. In fact, 
pQCD underestimates the p/tt + ratio by a factor of 3-10 
in pp collisions (see Fig. 3(a) of this work for s/s = 27.4 
GeV, and e.g. Ref. || for Tevatron energy). 

In this Letter we look into how pQCD is used to calcu- 
late px spectra in pp collisions. We focus on the role of a 
non-perturbative ingredient, the fragmentation function 
(obtained by fitting data) in the production of hadronic 
final states. Most of the information included in the fits 



comes from h + + h~ data. The fragmentation function 
(FF) of pions is studied in some detail. Less direct in- 
formation is available on kaons, and the FF of protons 
is even less well-known. In the following, we concentrate 
on the proton FF as an example of the role of the FF in 
the hadroproduction of hard particles. We find (for all 
types of hadrons) that the value of the FF in a region of 
phase space where it is least constrained plays a decisive 
role for hadron production at RHIC and LHC. 

To predict the px spectra of final state hadrons in the 
framework of pQCD, perturbative partonic cross sections 
need to be convoluted with parton distribution functions 
(PDF-s) and fragmentation functions (FF-s) according 
to the factorization theorem JjJ. Perturbative QCD has 
nothing to say about the details of FF-s, apart from de- 
scribing their scale evolution. The FF-s are assumed to 
be universal, meaning that once extracted from a limited 
set of data via a fitting procedure, they can then be used 
with predictive power in other reactions [p[-pd|. The sit- 
uation is conceptually identical to that of the PDF-s. We 
will work with this formalism in the present paper. As 
an alternative, phcnomcnological models (e.g. string or 
cluster models) can of course be implemented in Monte 
Carlo programs |l2|,f3) ■ 



The available fragmentation functions are obtained by 
fitting mostly e + e~ data. The FF-s extracted this way 
are preferred to using information also from pp (pp) data 
to avoid the complication introduced by kx smearing 
|p^ , ^5| in the latter case. Note, however, that the phase 
space for hadronic collisions is different from that of e + e~ 
reactions. This manifests itself in different values for z, 
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the momentum fraction of the fragmenting parton car- 
ried by the final hadron. In order to study the role of 
different regions of z in the FF, we define the ratio R z as 



R 2 (p T )= fdz' 

Jo prdprdz' 



, / da 
p T dp T dz' 



(1) 



where the cross section da/(prdpTdz) is differential in 
Pt and in z, and the ratio depends on the upper limit of 
the integration in the numerator. 



also very important for the hadroproduction of pions and 
kaons, which have harder FF-s. Let us examine, there- 
fore, what kind of constraints are placed on the FF-s in 
the large z region by e + e~ data. 

In e + e~ collisions, the hadron production cross section 
(differential in incoming momentum fraction x) is given 

by! 



da h 
dx 



Yl J -j-^-^y^F,^R)D h a {z^ F ) 



(2) 
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FIG. 1. The ratio Rz(pr) defined in Eq. (1) as a function 
of z for different pr values (a) at yfs = 27 A GeV with pr — 2 
GeV (solid line), p T = 4 GeV (dashed), and p T = 6 GeV 
(dotted), and (b) at y/s = 130 GeV with p T = 4 GeV (solid), 
p T = 8 GeV (dashed), and p T = 16 GeV (dotted). 

Figure 1 displays R z (pt) as a function of z for three dif- 
ferent values of pt in the case of proton production in pp 
collisions at center-of-mass (cm) energies y/s = 27 A GeV 
(part (a)) and y/s = 130 GeV (part (b)). In Fig. 1(a) 
we show pr — 2 GeV (solid line), 4 GeV (dashed) and 
6 GeV (dotted). In Fig. 1(b), where the higher cm en- 
ergy allows higher pt values, the solid line corresponds 
to pt = 4 GeV, dashed means pt = 8 GeV, and dotted 
represents pt = 16 GeV. It can be read off from Fig. 1(a) 
that when pt = 4 GeV, about 90% of the cross section 
comes from the contribution of z > 0.6 at y/s — 27 A 
GeV. When pt = 6 GeV, 90% of the cross section arises 
from z > 0.7. Fig. 1(b) shows that at y/s = 130 GeV, 
when pt = 4 GeV, approximately 60% of the cross sec- 
tion comes from the contribution of z > 0.5; when pt = 8 
(16) GeV, 70% (95%) of the cross section is from z > 0.5. 

In summary, the large z part of the fragmentation func- 
tion dominates the proton production cross section in 
hadronic collisions. When pt increases, the effect be- 
comes even more pronounced. This conclusion is consis- 
tent with results for Drell-Yan processes jl(| . 

The above results depend on the shape of the FF. If 
the corresponding FF-s fall less steeply (are "harder") 
than the KKP proton FF-s, the dominance of the large z 
region is even more significant. Thus, the above effect is 



where da a /dy is the differential cross section of a par- 
tonic sub-process producing parton a (as a function of 
y = x/z, of the factorization scale [ip, and of the renor- 
malization scale hr), and D^(z, /if) is the FF for parton 
a to fragment into hadron h with momentum fraction z. 
The contributions from different partons (quarks of flavor 
i and gluons g) are summed over. 

Through next-to-leading order (NLO), in the MS 
scheme, the partonic sub-process cross section for quark 
flavor i can be written as S 



da, 



dtj 



' -<r N c e qz 2 



6(1 -y) + 



2n 



(3) 



where ao is the corresponding total cross section for 
e + e~ — * 7* —>■ /i + /i~ , N c is the number of colors, e 9i 
is the charge of the quark, a s is the strong coupling con- 
stant, and the interested reader is referred to || for the 
functions P and K. The corresponding expression in- 
cluding the Z° contribution is similar, but more compli- 
cated due to the weak coupling of the Z° [Q . 
For gluons, 



P 



H'F 



+Kj(y)+K^y) 



(4) 



in a similar notation (Nf is the number of active flavors). 

The usual parameterization of the fragmentation func- 
tions at an input scale //q is p| [TT| 



D(z,nl) =Nz a {l-z) 







(5) 



where a and (3 are fixed by fitting to a given set of data. 
The parameter a describes the behavior of the FF in the 
small z region, while (3 determines the behavior in the 
large z region. The FF-s are not as well studied as the 
PDF-s. The KKP fragmentation functions || provide 
one of the few sets that contain FF-s for protons. Since 
we are particularly interested in proton production, we 
will use KKP FF-s || in the following calculations. Fig- 
ure 2 illustrates how arbitrary changes in the value of (3 
modify the proton FF-s in e + e~ collisions in the large z 
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region, and how these modifications influence the agree- 
ment with the experimental data at yfs = 91 GcV JL8| . 

In the leading order approximation, one has x = z. 
Therefore, at sufficiently high cm energies, e + e~ data 
points at large x should constrain the fragmentation func- 
tion at large z. It can be seen from Fig. 2 that the num- 
ber of available data points at large x is quite limited, 
and that the uncertainty associated with the highest-x 
point in particular is rather large. This is not surprising, 
since this point is close to the phase space limit of the ex- 
periment. Since there is no reason to expect that a NLO 
calculation would change the state of affairs at the high 
energies considered here, we use the leading order pQCD 
results in the followingdiscussion. We use [If — \fs (this 
sets the log terms in (g) and (^) to zero). 




FIG. 2. Perturbative QCD calculations using FF-s with dif- 
ferent large z behavior compared to e + e~ data at y/s = 91 
GeV JTI|. The solid line is the result with fi = (3 K kp [§. The 
dashed (dotted) line represents (3 = 0.8/3kkp (/? = 0.5/3xkp) 
for modified u, d (u, d) quark and gluon FF-s. 

Figure 2 illustrates that there is significant freedom in 
the choice of the large z behavior of FF-s based on e + e~ 
experiments. To clearly display the point that these data 
can be fitted by FF-s being very different in the large 
z region, we show the cross section calculated with the 
KKP value § of /3 in Eq. (§) (solid line), together with 
results with (3 = 0.8(3kkp (dashed) and (3 = 0.5Pkkp- 
We conclude from Fig. 2 that taking (3 = 0.&(3kkp does 
not change the quality of the fit to the e + e~ data at this 
energy. It is also clear from the above discussion, how 



different parameterizations of FF-s |Sj 1 1 can be rather 
different in the large z region even for hr + h~ fragmen- 
tation (the FF-s studied in most detail). 

The situation concerning gluon FF-s at large z is even 
less certain than it is for quarks. The contribution of 
gluon fragmentation to the e + e~ hadron production cross 
section appears only as a NLO correction. In hadronic 
collision, the contribution from gluon processes is much 
more important than in e + e~ Jl6[ . Since the probability 
of finding a gluon in the proton (the gluon PDF) increases 



rapidly as x ~ 2pr/\/s decreases, gluon fragmentation 
plays an amplified role at RHIC and LHC energies, com- 
pared to fixed target energies. It is thus very important 
to note that the large z gluon FF for h + + hr obtained 
including some pp data can be an order of magnitude 
larger than the one not including pp information jllj . 

As we have seen, the e + e~ data do not strongly con- 
strain the FF-s in the large z region. It is therefore felt 
that one has some freedom to fit the proton px spectra in 
pp collisions by modifying the large z behavior of the pro- 
ton fragmentation functions. To illustrate this idea, we 
compare the results for the p/tt + ratio calculated using 
the KKP form of the proton FF-s, but varying the value 
of the parameter (3. In addition, the transverse momen- 
tum distribution of the partons in the proton needs to be 
taken into account in a more complete calculation. As 
long as no first-principles treatment of this effect is avail- 
able, the width of the transverse-momentum distribution 
(fef.) provides another phenomenologically adjustable pa- 
rameter. 




7 8 
p T (GeV) 

FIG. 3. Comparison oip/n ratio data and pQCD results in 
pp collisions: (a) at yfs — 27 .4 GeV ( dotted line is j3 = [3k kp, 
dashed is j3 = 0.8/3kkp, solid line is calculated with different 
(fcy) values for the pion and the proton as described in the 
text; data are from Ref. Q), (b) at y/s = 61 GeV (ISR data 
§2)) and prediction for RHIC at yfs = 130 GeV (dashed line). 

Figure 3 demonstrates that the remaining pQCD free- 
dom represented by the above parameters may be suffi- 
cient to achieve a better description of the available data 
on p/7r + ratios, without invoking any other mechanism. 
Our goal here is not to fit the data; Fig. 3 serves only as 
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an illustration. The dotted line in Fig. 3(a) shows that 
the calculated R — p/ir + ratio using the original KKP 
proton FF-s underestimates the data (dots) at y/s = 
27 A GeV by up to a factor of 10. However, if we set 
= 0.80k kp in the proton FF (dashed line), the pQCD 
result for the p/n + ratio comes close to the experimental 
data for px > 6 GeV. In the above calculation, we used 
(k T ) = 0.9 GeV 2 for both proton and pion production. 
However, it is not clear that the input (fc 2 ^) for proton and 
pion production needs to be identical. Due to the large 
mass of the proton, the effective s (the energy involved 
in the partonic cross section) is larger for proton produc- 
tion than for pion production at the same px- Large s 
leads to more room for a dynamical intrinsic fc-r, just like 
in the Drell-Yan case, where the larger Q 2 of the lepton 
pair leads to a larger k T . Therefore it is conceivable that 
a larger input (k 2 -) is needed for proton production than 
for pion production p0| . 

Using (k T ) p = 2.6 GeV 2 for the proton and keeping the 
value (fcf^Tr = 0.9 GeV 2 for the pion in the calculation, 
we obtain the solid line in Fig. 3(a) for the p/ir + ratio 
in pp collisions. It is fair to say that the result of this 
calculation is in satisfactory agreement with the data for 
p T > 3 GeV. (Below p T « 2-3 GeV pQCD is not consid- 
ered reliable, and non-perturbative effects may become 
important pl|.) We obtain similar results comparing to 
data at = 38.8 GeV. 

In Fig. 3(b) we explore the energy dependence of the 
above proposition. Here, the solid line is the result of a 
pQCD calculation for the p/ir + ratio in pp collisions at 
y/s = 61 GeV, with = O.80 KKP and (k^) p = 3 GeV 2 
(while keeping (k T ) n constant). The calculated results 
are compared to ISR data (open symbols) [g2j. We also 
give the prediction for RHIC energy, y/s = 130 GeV, with 
(fc 2 ) p = 3 GeV 2 (dashed line). 

The main purpose of this paper is to call attention to 
the importance of the behavior of fragmentation func- 
tions in the large z region for hard hadron produc- 
tion. Recently, similar ideas were considered for bot- 
tom hadroproduction at the Tevatron in terms of the 
moments of the fragmentation function (higher moments 
emphasizing the importance of large z) p3fl . 

In conclusion, since the large z details of the fragmen- 
tation functions are very important for pQCD predictions 
of high pt hadron production in both pp and AA colli- 
sions at RHIC and LHC, a comprehensive study of these 
fragmentation functions with particular attention to the 
large z region is strongly warranted. Additional com- 
plications not addressed in the present note include e.g. 
the effect of changing on the scale dependence, a first- 
principles treatment of the transverse momentum degree 
of freedom, and various nuclear effects. These are left for 



future work. It appears that available p/ir ratios from pp 
collisions at high px can be reproduced by adjusting the 
large z behavior of the fragmentation functions and the 
width of transverse momentum distributions. 
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